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ABSTRACT 

We show that the winds identified with high-redshift low-mass galaxies may strongly affect the for- 
mation of stars in more massive galaxies that form later. With 3D realizations of a simple nonlinear 
growth model we track gas shocking, metal enrichment, and cooling, together with dark halo formation. 
We show that outflows typically strip baryonic material out of pre-virialized intermediate mass halos, 
suppressing star formation. More massive halos can trap the heated gas but collapse later, leading to 
a broad bimodal rcdshift distribution, with a larger characteristic mass and metallicity associated with 
the lower redshift peak. This scenario accounts for the observed bell-shaped luminosity function of 
early-type galaxies, explains the small number of Milky- Way satellite galaxies relative to standard Cold 
Dark Matter proscriptions, and provides a reasonable explanation for the lack of metal poor stars in the 
solar neighborhood and the more general lack of low-metallicity stars in other massive galaxies relative 
to "closed-box" models of chemical enrichment. Heating of the intergalactic medium by early outflows 
should produce spectral distortions in the cosmic microwave background that will be measurable with 
the next generation of experiments. 

Subject headings: galaxies: formation - galaxies: interactions- large-scale structure of the universe - 
galaxies: elliptical 



1. INTRODUCTION 

It has long been recognized that the X-ray luminosity- 
temperature (L X ~T) relation of clusters does not obey the 
simple scaling laws that would hold if clusters were formed 
from the collapse of unheated primordial gas, and thus the 
gas within clusters is likely to have been heated before the 
formation of the clusters themselves (Kaiser 1991). Sub- 
sequent investigations have determined that preheating is 
also necessary to explain the L x — T relation as a function 
of cluster mass (Cavaliere, Menci, & Tozzi 1999). 

The high metallicity of cluster gas and the claimed over- 
abundance of alpha elements (Gibson and Matteucci 1996; 
Lowenstein & Mushotzky 1996) points to preheating by 
SNII driven winds (Renzini et al. 1993; Trentham 1994; 
Nath & Chiba 1995). High enrichment is not restricted to 
the most massive clusters, but appears to be widespread, 
extending to groups of galaxies (Buote 2000). Empirically 
the epoch of this enrichment is now known to occur at 
z > 0.5, and approximately solar enriched cluster gas has 
been detected at redshifts as large as z ~ 1 (Hattori et al. 
1997). Hence it is natural to suppose that pre-enrichment 
and preheating are the consequence of very early and vig- 
orous massive star-formation. 

Cluster gas is evidently a sink for enriched and heated 
material, and the most likely culprit for this enrichment 
is dwarf galaxies. Theoretical work has shown that super- 
novae and OB winds in these low-mass objects should lead 
to the production of energetic outflows with temperatures 
on the order of 10 6 K (Larson 1974; Dekel & Silk 1986; 
Vader 1986). This behavior has been clearly identified in 
studies of both local starbursting galaxies (Axon & Taylor 
1978; Marlowe et al. 1995; Heckman 1997; Hunter et al. 
1998; Martin 1998) and spectroscopy of high-z galaxies ( 



Franx et al. 1997; Pettini et al. 1998; Frye & Broadhurst 
1998; Warren et al. 1998). Whether these outflows lead to 
a catastrophic loss of the interstellar gas however, is likely 
to depend on a number of factors (De Young & Heckman 
1994), and is a subject of current investigations (Murakami 
& Babul 1999; Mac Low & Ferrara 1999; Strickland & 
Stevens 1999). 

In the generally investigated hierarchical models of 
structure formation such as the Cold Dark Matter model 
(CDM), the existence of an era of widespread enrichment 
by outflows from dwarf galaxies is in fact quite natural, 
because low-mass galaxies are expected to form in large 
numbers and at early times (e.g., White & Frenk 1991). 
Since these early galaxies will be found preferentially in 
the large-scale overdense regions that later form clusters, 
dwarf outflows are the obvious candidates for pre-cluster 
heating. 

The existence of large numbers of small galaxies at 
high redshifts is also favored by observations, to help un- 
derstand the steep number counts and low luminosities 
of faint galaxies (Broadhurst, Ellis, & Glazebrook 1992), 
the sizes of faint galaxies in Hubble Deep Field images 
(Bouwens, Broadhurst, & Silk 1998a, b), and the small 
sizes of the distant Lyman-break galaxies (Steidel et al. 
1999). Locally, however, the space density of dwarf galax- 
ies relative to massive galaxies is far less than predicted 
on the basis of the steep Press-Schechter slope for the 
faint end of the mass function (Ferguson & Binggeli 1994), 
prompting theoretical studies of the disruption of dwarves 
by tidal forces from neighboring objects (Moore et al. 
1998), external UV radiation (Kepner, Babul, & Spergel 
1997; Norman & Spaans 1997; Corbelli, Galli, & Palla 
1997), and catastrophic mass loss during outflows (Larson 
1974; Dekel & Silk 1986; Vader 1986). 
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While many of these mechanisms may have had an im- 
pact on the formation of dwarf galaxies, relatively lit- 
tle attention has been directed towards the influence of 
outflows on neighboring galaxies. As the earliest galax- 
ies to form were highly clustered (Kaiser 1984) and typi- 
cal outflow temperatures and velocities were much larger 
than the virial temperatures and velocities of these galax- 
ies, it is likely that dwarf galaxies were strongly influ- 
enced by their neighbors (Scannapieco, Ferrara, & Broad- 
hurst 2000). Similarly, the sources responsible for pre- 
enrichment of the intracluster medium may well have en- 
riched larger protogalaxies, with important consequences 
for their metallicity histories and cooling times. 

In this work, we conduct an idealized investigation of the 
impact of dwarf outflows on the history of galaxy forma- 
tion. While the consequences of homogeneous heating on 
galaxy formation have been examined in the past (Blan- 
chard, Valls-Gabaud, & Mamon 1992), the inhomogeneous 
nature of this process and the associated pre-enrichment 
of galaxies have not been addressed. Our aim is not, how- 
ever, to construct a complete model of galaxy formation, 
and thus we do not track processes such as the forma- 
tion of second-generation stars in galaxies, production of 
dust, transfer of angular momentum, or the structure of 
the interstellar medium. Rather we focus on the properties 
of the inhomogeneously heated and enriched intergalactic 
medium (IGM) out of which galaxies coalesced, the likely 
consequences of galaxy formation in this environment, and 
to what degree these issues must be accounted for within 
more detailed simulations of galaxy formation. 

The structure of this work is as follows. In §2 we re- 
view the observational evidence that leads us to consider 
a model in which widespread dwarf outflows shocked and 
enriched the medium out of which larger galaxies formed. 
In §3 we describe a simple numerical Monte Carlo code 
that we use to asses the overall features of such a model. 
The results of our simulations are given in §4 in which 
we examine which aspects of galaxy formation are most 
sensitive to the presence of outflows. In §5 we discuss the 
limitations of our modeling, and conclusions are listed in 
§6- 

2. OBSERVATIONAL EVIDENCE 

In this section we review the observational evidence 
that points towards widespread early enrichment by dwarf 
galaxies. We show that the chemical and thermal proper- 
ties of nearby galaxy clusters are both suggestive of such 
pre-enrichment, as is the weak evolution of these properties 
with redshift. We show that such a picture is consistent 
with our knowledge of the intergalactic medium as well, 
and helps to explain the relatively high metallicity seen 
in many Ly absorption systems. The presence of outflows 
has been identified in emission line studies of high-z galax- 
ies, and several nearby dwarf galaxies have been caught 
"red-handed" , surrounded by clouds of enriched gas with 
temperatures greatly exceeding virial. Finally, we review 
some unsolved questions in galaxy formation that may de- 
pend on preheating and enrichment and serve to motivate 
our exploratory numerical studies. 

2.1. Galactic Outflows and the Intra Cluster Medium 

Many cluster properties can be understood in the con- 
text of self-similar models (Kaiser 1986, 1990) which are 



exact for power-law initial fluctuation spectra and a rea- 
sonable approximation for more realistic spectra such as 
that of the Cold Dark Matter (CDM) model. While the 
optical properties of galaxies in clusters are consistent 
with these predictions, the slope of the X-ray luminosity- 
temperature relationship of the hot intra-cluster gas is too 
steep and slowly evolving to be understood in this context. 

A convincing explanation of L x — T evolution has been 
made by Kaiser (1991) who speculated that at an epoch 
before the formation of present-day clusters, the gas which 
formed the intracluster medium (ICM) was preheated, in- 
jecting sufficient energy to expel gas from the small po- 
tential wells that were non-linear at early times. By as- 
suming that this heated gas later cooled adiabatically onto 
the large cluster-size potential wells, he was able to sim- 
ply reproduce the observed X-ray properties of clusters. 
Numerous investigations of preheating and L x — T evolu- 
tion have also reached similar conclusions (David, Jones, 
& Forman, 1996; Mushotzky & Scharf 1997; Eke, Navarro, 
Frenk 1998) as have studies of the L x — T relation in groups 
of galaxies (Cavaliere, Menci, & Tozzi 1999). 

The origin of this preheating is most naturally due to 
Supernovae Type-II (SNell) activity in the dwarf star- 
burst population. Several studies of expanding HI gas in 
nearby dwarf galaxies show clear evidence of dense ex- 
panding shells with velocities above 15 km/s (Marlowe et 
al. 1995; Heckman 1997; Hunter et al. 1998; Martin 1998). 
Lya studies confirm that similar outflows are present in 
higher z galaxies at rcdshifts of order <~ 3 (Pettini et al. 
1998) and even higher (Frye & Broadhurst 1998; Warren 
et al. 1998 ), where the overall space density of small star- 
bursting galaxies is larger. 

A natural prediction of such a picture would not only 
be the heating of the ICM by the dwarf population, but 
the enrichment of this gas by metals expelled by the SNell 
powering the galactic winds. In fact, the metallicity of the 
intracluster medium is observed to be quite high and con- 
stant (w 0.3Zq) over a large range of cluster masses (Ren- 
zini 1997). This roughly constant value is evidence that 
the clusters have neither lost or gained a large amount of 
baryons during their evolution, as massive cluster outflows 
or accretion of pristine gas would lead to a large scatter. 
Furthermore, this metallicity is roughly the same today as 
observed at z w 0.3 and perhaps even at z w 1 (Hattori 
et al. 1997), suggesting that IGM enrichment took place 
early on in the lifetime of the clusters (Renzini 1999). 

Several authors have investigated the possibility that 
this enrichment is due to the present cluster members and 
that most of the X-ray emitting hot gas in clusters is due to 
outflows from galaxies that populate them today. Okazaki 
et al. (1993) considered gas ejection from bright, elliptical 
galaxies in clusters, concluding that these galaxies could 
contribute no more than 10% of the ICM. Trentham (1994) 
suggested that the ICM instead may be primarily due out- 
flows from precursors to the observed populations of dwarf 
galaxies in clusters, but Nath and Chiba (1995) concluded 
that the metallicities generated in this type of scenario 
were sufficient to explain only clusters with low-metallicity 
gas. 

Perhaps most convincingly Gibson and Matteucci (1996) 
have studied the possibility that the elliptical galaxies ob- 
served in clusters today were responsible for the major- 
ity of ICM enrichment. By developing models of galac- 
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tic winds consistent with the observed properties of clus- 
ter ellipticals they concluded that even in their "maximal 
models" in which all the gas returned by dying stars is 
ejected into the ICM, neither the giant elliptical, the dwarf 
spheroidals, or both these populations combined can be re- 
sponsible for more that 33/38% of the ICM. And yet this 
gas is hot and metal rich. 

A final piece of evidence as to the origin of ICM metals 
comes from studies of relative abundances. While for many 
years, studies of metals in clusters were limited to measure- 
ments of iron abundances, the X-ray Advanced Satellite for 
Cosmology and Astrophysics (ASCA) (Tanaka, Inoue,& 
Holt 1994) provided the first opportunity to study other 
metals. This allowed Loewenstein & Mushotzky (1996) 
to observe excesses of alpha-elements such as O, Mg, and 
Si, which indicate that most of the ICM enrichment was 
due to SNell rather than Type la (SNela). This result 
was later contested by Ishimaru & Arimoto (1997), who 
claimed that SNela enrichment could be responsible for 
50% or more of ICM metals, based on the solar 'mete- 
oritic' metallicity rather than the solar coronal gas. Gib- 
son, Loewenstein, & Mushotzky (1997) pointed out, how- 
ever, that this result was specifically linked to the SNell 
yield used in their analysis. While they were unable to rule 
out Ishimaru & Arimoto's claim within theoretical uncer- 
tainties, they showed that more recent SN models that 
treat convection and mass-loss, reduce the fractional con- 
tribution ICM by SNela to less than 5%. Renzini (1997) 
on the other hand while confirming the importance of early 
Snell, has argued that later enrichment by significant num- 
bers of SNela are required to produce the observed high 
level of Fe relative to alpha elements. 

In summary then, metallicities studies of clusters indi- 
cate that the ICM was enriched early on, that this enrich- 
ment is too high to be due to purely to the observed cluster 
members, and that the relative abundances of metals are 
suggestive of widespread SNell enrichment; all indicative 
of enrichment by high-redshift dwarf galaxies. 

2.2. Galactic Outflows and the Intergalactic Medium 

The first galaxies to form would have necessarily had a 
huge impact not only in regions that would later form clus- 
ters, but on the Intergalactic Medium (ICM) as a whole. 
In standard schemes for hierarchical structure formation, 
the first baryonic objects formed at redshifts ~ 40 and with 
masses of order 1O 5 M , somewhat smaller than the dwarf 
galaxies observed today (Haiman, Thoul, & Loeb 1996). 
While several authors have shown that the total flux of 
UV photons from stars within these objects is sufficient to 
reionize the universe (Couchman & Rees 1986; Fukugita 
& Kawakasi 1994; Shapiro, Giroux, & Babul 1994), these 
objects are likely to suppress their own formation long 
before this occurs. As molecular hydrogen is easily photo- 
dissociated by 11.2-13.6 eV photons, to which the uni- 
verse is otherwise transparent, the emission from the first 
stars quickly destroys all avenues for cooling by molecu- 
lar line emission. This quickly raises the minimum virial 
temperature necessary to cool effectively to approximately 
10,000 K, suppressing the further formation of objects with 
masses <; 10 8 M© (Haiman Rees, & Loeb 1997; Ciardi et 
al. 1999). 

Thus reionization, if achieved by galaxies, was relatively 
late (z < 15) and associated with objects of similar size 



or larger than dwarf galaxies. This high mass scale, im- 
plies that the IGM phase transition must have occurred 
in a clumpy and inhomogeneous manner, with impor- 
tant implications for small-scale microwave background 
anisotropies (see eg. Aghanim et al. 1996; Miralda-Escude, 
Haehnelt, & Rees 2000; Scannapieco 2000). 

Early enrichment and heating is similarly expected to 
be inhomogeneous, following the spatial distribution of 
the first galaxies which are restricted to rare overdense 
regions. Hence measurements of a low average metallicity 
(~ O.OIZq) of the Ly-forest at high redshift (Songalia & 
Cowie 1996; Savaglio 1997; Songalia 1997) may be viewed 
not so much as evidence for a low metallicity intergalactic 
medium, as an inhomogeneous one. High-redshift regions 
of higher density that later develop into clusters would 
therefore represent the most polluted volumes, where one 
would naturally expect a significantly higher mean metal- 
licity. 

Direct evidence that the IGM is enriched by galaxy out- 
flows may be inferred from detailed observations of local 
dwarf spheroidal galaxies. For example, ASCA observa- 
tions (della Ceca et al. 1996) show that the star-forming 
dwarf galaxy NGC 1569 is surrounded by a hot (8 x 10 6 K) 
halo of gas whose temperature greatly exceeds the virial 
temperature of the galaxy and whose line strengths are 
consistent with O.25Z0. Similarly, X-ray observations of 
the dwarf irregular NGC 4449 by Bomans, Chu, & Hopp 
(1997) show that it is embedded in a supergiant shell of 
~2x 10 6 K gas with a metallicity of ~ 0.3Zq. 

2.3. Galactic Outflows and Galaxy Formation 

The presence of an inhomogeneous IGM with strong 
temperature and metallicity fluctuations creates quite a 
different environment for galaxy formation than the pri- 
mordial conditions often assumed. This becomes clear 
when one contrasts the virialization and cooling of pri- 
mordial gas to that of gas pre-enriched and heated by 
dwarf outflows. In Figure |l| we replot a classic compari- 
son first made in Dekel and Silk (1986). Collapse redshifts 
of lcr to 2cr dark matter halos (dotted lines) are calcu- 
lated from linear theory for both a flat CDM (f2o = 1, 
n A = 0,n b = 0.07, cr 8 = 0.6, and T = 0.44) and a ACDM 
model (fi = 0.35, Q A = 0.65, Q b = 0.06, <j 8 = 0.87 and 
r = 0.18), assuming spherical collapse. The solid lines 
show the lowest redshift that a sphere of gas can virialize 
and still have time to cool and form a galaxy by z = 0, 
as calculated by the cooling models described in §3.2. In 
each pair the upper line corresponds to primordial gas, 
and the lower line to gas that has been enriched to OAZq. 
Here we see that even modest enrichment by outflows has 
the potential to greatly accelerate galaxy formation on the 
~ 10 12 M & /h scale due to the long cooling times of large 
clouds without the additional avenue for cooling afforded 
by line emission by metals. 

Also on this plot, we show the mass of halos associated 
with a virial temperature of 5 x 10 5 K, typical of galaxy 
outflows (dashed line). This serves as an estimate of the 
smallest mass of galaxies that can form in areas impacted 
by outflows without being disrupted by shocks. This sug- 
gests that while high-mass galaxy formation may be en- 
hanced by dwarf outflows, the formation of dwarf galaxies 
has the potential to be suppressed. 
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The features shown on this plot invite comparison with 
the properties observed in elliptical galaxies, which are 
preferentially found in the most enriched regions in the 
universe, and whose mass functions are biased to large 
values relative to the rest of the galaxy population. It is 
also interesting to note that various analytical and N-body 
studies of cold dark matter models (Kauffmann, White, & 
Guiderdoni 1993; Klypin et al. 1999; Moore et al. 1999) 
have shown that ~ 50 satellites with circular velocities 
~ 20 km/s should be found within 600 kpc of the Galaxy, 
while only 11 are observed. Again, this lack of local dwarf 
galaxies is suggestive of galaxy formation in an intergalac- 
tic medium that has been impacted by outflows. 

Closer to home, the abundance distribution of long- 
lived stars in the solar neighborhood shows far too few 
low-metallicity stars compared with a simple "closed box" 
model of galactic chemical evolution. This is the long 
standing G-dwarf problem, first pointed out by van de 
Bergh (1962) and Schmidt (1963). The sudden drop in 
the number of G-dwarfs with metallicities below ~ O.IZq 
has lead to a number of models in which an initial pro- 
duction spike or "prompt initial enrichment" adds metals 
to the gas out of which the majority of stars form (e.g., 
Truran & Cameron 1971; Ostriker & Thuan 1975: Koppen 
& Arimoto 1990). This ad hoc floor to the metallicity then 
allows good fits to the Galactic stellar data. 

Similar pre-enrichment may also be necessary to explain 
the metallicity distribution of stars outside our own galaxy. 
Integrated spectra of elliptical galaxies have been studied 
with care (Worthey, Dorman, & Jones 1996) concluding 
rather puzzlingly that such galaxies do not have a closed- 
box history but again show the need for a "floor" of around 
0.1Z e (Thomas, Greggio, & Bender 1999). 

3. COSMOLOGICAL MONTE CARLO 

The great preponderance of observational clues and 
pointers suggests that we should seriously consider a sce- 
nario in which the formation of modern-day galaxies was 
preceded by an era of IGM enrichment by a high-redshift 
population of dwarf starbursting galaxies. Note that from 
a theoretical point a view, such a model is not ad hoc, but 
rather a consequence of hicrarchal structure formation. As 
smaller objects form early and are able to suppresses fur- 
ther small-scale formation, it is only natural that larger 
galaxies would have formed in a second wave of collapse 
within this pre-enriched medium. 

The density dependence of the effect of outflows on 
galaxy formation requires a 3-D calculation. To date 
only the average properties of the IGM as a function of 
time have been explored in one dimensional calculations 
by Blanchard, Valls-Gabaud, & Mamon (1992) and by 
Nath and Trentham (1997). Such averages are appropri- 
ate to the IGM as a whole, but underestimate the effect 
on galaxy formation. Thus for example, the outflow model 
proposed by Nath and Trentham (1997) produces metal- 
licities ~ .01Z & , while the observed metallicities of cluster 
gas and elliptical galaxies are comparable to solar abun- 
dance (Renzini 1993). 

In order to determine the impact of outflows on galaxy 
formation, we have developed a simple cosmological Monte 
Carlo code, in which we realize a small volume of the uni- 
verse and study the linear evolution of objects in the mass 
range 2 x 10 8 n /hM Q < M < 1 x 10 13 n /hM Q within it, 



with reasonable extensions to deal with collapse. Our phi- 
losophy is not to pretend that we can possibly reproduce 
in detail the complicated processes of dark matter halo col- 
lapse, galaxy formation, and gas infall and expulsion, but 
rather to explore a model that captures the essential fea- 
tures and predicts with some confidence rough magnitudes 
and trends that can be compared with observations. 

3.1. Collapse of Dark Matter Halos 

We choose for our simulation a cubic comoving volume 
of (12Mpc//i -1 ) 3 divided into 512 3 cells and with periodic 
boundary conditions. On this mesh we construct a cos- 
mological linear overdensity field <5(x, z) — p(x, z)j p(z) 
where p(z) is global average density of the universe, x is 
a position in comoving coordinates, and z is the redshift. 
Transforming into Fourier space the density field can be re- 
expressed as <5(k, z) = J (i 3 xexp(— ik-x) <5(x, z). From lin- 
ear theory, the evolution of the density field as a function 
of time is given simply by <5(k, z) = 8qQ&)D{z) / Dq where 
D(z) is the dimensionless growth factor, S (k) = <5(k, 0), 
and D = D(0). 

In the CDM model, the Gaussian random field S(k) can 
by constructed by randomly choosing modes such that 
the variance is given by (<5(k)<5(k')) = (5(k)5*(— k')) = 
(27r) 3 <5(k + k')P(k). The power spectrum is then 



P(k) = 2Tr 2 Sj { k n T 2 (k p Mpc/hT), 



(1) 



where T(q) is the CDM transfer function given by Eq. 
(G3) of Bardeen et al. (1996), and we choose aoH = c, 
such that k p = k/a = kH /c is the physical wave 
number. The shape parameter T is defined as T = 
O /icxp(— Clb — ^b/^o) (Sugiyama 1995) and is observed 
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given a spectral index of n 



(Vianna & Liddle 1996). We fix the normalization fac- 
tor, 5h, against the number abundance of clusters. The 
variance of the mass inclosed in a sphere of radius R is 
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If we choose a spherical top-hat window function defined 

then we can normalize our 
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by W{x) = 3 
fluctuation spectrum by setting a% = <r(8Mpc/r 
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in the open case and C(f2o) = -59 + O.16f2o - 
take (fi + = 1) (Vianna & Liddle 1996) 

Having constructed <5o(k) in this manner, we then 
convolve it with window functions with lengths scales 
corresponding to ten masses arranged logarithmi- 
cally from Mi = 2 x 10 8 O /hM Q to M w = 
10 13 ^ / 7iM , to obtain ^(x) = 5 (x,R2xio»n /hM ) = 
I 757^ exp(ik • x) S(k)W(R 2xW *n /hM @ k) 



<5,f 2 (x) 



5 (x, R eS xwHi„/hM e ), ■■■ ^"W. tnc linear density field 
at z — smoothed at each mass scale. Here the mini- 
mum mass of our simulations is set by molecular cooling 
constraints discussed further below. 

From the spherical collapse model we can identify the 
linear overdensity <5 sc (z) at which the true density field 
has virializcd. For flat models this value is ~ 1.69 at all 
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times while in open cases it is a weakly decreasing function 
of redshift as fitted in Appendix A of Kitayama & Suto 
(1996). For each peak in the linear overdensity held such 



that S i (x, R) > S sc (z = 0), we identify a collapse redshift 
z= r such that 



£jj**(x) = S c (z sc ) 



D 



D{z s , 



(3) 



This approach has the limitation that it becomes inac- 
curate as <r(R)D(z)/D approaches <5 SC . Thus while the 
spherical collapse model works well at determining the 
number and spatial clustering of rare objects (Lacy & Cole 
1993; Mo, Jing, & White 1996, 1997) it fails to agree with 
numerical simulations for more common low-mass halos 
(Lacey & Cole 1994; Sheth & Tormen 1999). Further- 
more, the object-by-object identification of linear peaks 
with collapsed objects has been shown to be unreliable, 
even in cases in which the statistical properties are in good 
agreement with N-body simulations (Bond et al. 1991). 

Sheth, Mo, & Tormen (1999) have shown that these pre- 
dictions can be improved by accounting for the ellipticity 
of collapsing clouds caused by tidal forces from nearby col- 
lapsing peaks. By including a correction factor for the crit- 
ical linear overdensity, they have been able not only to im- 
prove statistical predictions over the spherical model, but 
do well on an object-by-object basis. As our Monte-Carlo 
approach depends on this object-by-object identification, 
we can greatly improve the accuracy of our simulations by 
comparing each peak with S ec (z), such that 
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where (3 = 0.47 and 7 = 0.615. 

Having identified the collapsed peaks by either al- 
gorithm, we then arrange them in order of decreas- 
ing collapse redshift, z c = z sc or z cc , to obtain a 
list of candidate points and redshifts for collapsed ha- 
los at each of the different mass scales: Candidate = 

r Mi Mi 1 , z Mi = r Mi Mi i p 

t x l > x 2 !— / dnu ^"c, candidate — \ Z c, 1 > z c,2 ) ■ • ■ J • fIom 

this list we exclude unphysical points corresponding to all 
halos that collapse within a previously collapsed halo at 
the same or greater mass scale. That is we remove all 
points such that 



Mi < Mj, ||x 



M; 
k 



Mi 



< R Mj , and z^ < z^l , 



(5) 



where the distance is calculated accounting for the peri- 
odicity of the simulation volume. 

With these points excluded, we then have the full his- 
tory of the collapsing dark matter halos at each of the 
various mass scales: X Mi = {x.^ ,x. Mi , ...}, and Z^ li = 

{z^{, z^2i ■••}• The merger history of these objects is also 
simply calculated by searching for all points such that 



M L <M 3 , ||xf* 



-4^\\<R Mj , andz^ >z$. (6) 



Each such case indicates that the halo centered at x b 4 was 
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absorbed into the larger halo centered at x ; 3 at a redshift 

c M i 
0iz c,l- 



3.2. Comparison with Analytical Results 

At each redshift, the number density of collapsed halos 
in the simulation can be compared to analytical predic- 
tions. For the spherical collapse model this gives (Press & 
Schechter 1974) 



dn sc (M,z) 
dM 



exp 



2p{z) S sc (z)D 
M a(M) 2 D{z) 
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2a 2 (M)D(z) 2 J dM 



(7) 



Strictly speaking, Eq. (J7]) corresponds to the number of 
peaks as determined by a sharp fc-spacc filter, rather than 
the number determined by the top-hot filter used in our 
code. Equating these quantities has become somewhat of 
a common practice however, and has been shown to be a 
good estimate. 

In Figure ^ we compare the number density of objects 
in our simulation as function of redshift with the num- 
ber predicted by Eq. (0). As both the Press-Schechter 
integral, and our peak finding algorithm become inaccu- 
rate at late times, we exclude all peaks that collapse with 
a redshift lower that a minimum value, z m [ n (Mi) such 
that a(M l )D(z min (M i )) / D = 1.3S SC . We consider the 
ACDM model shown in Figure 0, (O = 0.35, ft A = 0.65, 
n b = 0.06, cr 8 = 0.87 and T = 0.18, h = 0.65). 

Here we see that the number of identified objects agrees 
well at most mass scales and redshifts. At the smaller 
scales, there is some discrepancy in the collapse time of 
objects, with the code slightly over-predicting the num- 
ber of objects formed at very early times. At larger-range 
scales, the Monte-Carlo approach does the best, matching 
the Press-Schechter predictions at all redshifts. Finally, at 
the highest mass scales, the low number of objects intro- 
duces significant statistical noise in our comparison. 

The spatial distribution of collapsed halos can also 
be compared to analytical predictions. If we define 
the correlation function at a smoothing scale R as 
the Fourier transform of the power spectrum smoothed 
by an appropriate top- hat window function, = 
2^ f °° k 2 dkP(k)W 2 {kR) s -^^-, then the correlation of 
collapsed peaks of scale R can be approximated by (Kaiser 
1984) 



1/2 



rfc(^/2 1 /2)] 2 x 



dye y2 / 2 erfc 



y-yU(r)/U(Q) 



y/2-2e R (r)/e R (0) 



(8) 



where erfc is the complimentary error function and v = 
5 c Do/cr(R)D(z). Note that this expression arises from con- 
sidering the fraction of Gaussian-distributed objects above 
a certain threshold, rather than from the more sophisti- 
cated excursion set approach (Bond, Efstathiou, & Kaiser 
1991; Lacy & Cole 1993) used to derive Eq. (§). 

In Figure || we compare the correlation function as given 
in Eq. (||) to that derived from the distances between the 
collapsed peaks at the five lowest mass-scales in our sim- 
ulation, with redshifts chosen such that v = 1.5 at each 
length scale. Eq. (j|) is intended to represent the corre- 
lations between collapsed objects at or above a particular 
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mass scale corresponding to the smoothing length scale. 
In practice, however, it makes little difference whether we 
compute £r m . (r) considering distances between all objects 
with masses above Mi or only between those in the mass 
bin corresponding to Mj itself. As comparisons only within 
a particular bin are both easier to calculate and under- 
stand physically, however, it is this quantity that we plot 
in Figure [|. 

Here we see that the numerical and analytical expres- 
sions agree at lengths scales ;> R Mi , but diverge at smaller 
distances. As the simulation simply excludes all points 
closer than Rm ( there is a sharp fall of in £r m (r) below 
this distance. The analytical expression, however, does not 
exclude pairs of objects that are contained within a larger 
virialized object (as would be excluded in an excursion set 
calculation) and increases dramatically at small r. Thus 
these quantities are derived in a manner such that they 
must be discrepant at small distances, and we can have 
some confidence that the overall distribution of objects in 
our simulations is reasonable. 

While the spherical collapse model is a logical testing 
ground for our code and can be simply compared to analyt- 
ical predictions, a more accurate simulation results by ap- 
plying the elliptical collapse criteria for halo formation, Eq. 
(Q) . This allows not only for a more accurate identification 
of collapsed peaks on an object-by-object basis, but delays 
the collapse of the more common peaks, and thus the min- 
imum accurate redshift at each mass scale is decreased to 
(Mi) such that a(M i )D(z ni i n (M i )) / D = 1.3<5 CC . 

In this model, the Press Schechter relation, Eq. (0), is 
modified to 



dn cc (M, z) 
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(9) 



where A = 0.322 and q = 0.3. Note that this function is 
essentially equivalent to the spherical collapse predictions 
for rare, high-mass peaks, and thus the ag normalization 
from the number density of galaxy clusters is unchanged 
in this model. This expression has been checked against 
the number density of objects in our simulation and again 
is in good agreement. 

3.3. Initial Gas Infall 

While our simple nonlinear collapse model captures the 
overall history of the dark matter distribution, the evolu- 
tion of the gas is more involved. While gas traces dark 
matter at early times, the formation of bound objects and 
subsequent star formation requires that the gas be able to 
cool to T ~ 0. 

If gas collapses and virializes along with a dark matter 
perturbation, and we assume an isothermal distribution, 
then it will be heated to a temperature of (Eke, Cole, & 
Frenk 1996) 



T ■ — 

-L vir — 



90K / 6.8 



f3 \5X + 3 



Af 6 2/3 (1 + z c ) 



tt Q A c (z c ) 
n(z c )18ir 2 
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( 10 ) 

where (3 is the ratio of specific galaxy kinetic energy to 
specific gas thermal energy, X is the hydrogen mass frac- 
tion which we take to be 0.76, M 6 = M/(l0 6 M Q /h), and 



A c is the ratio of the mean halo density to the critical 
density at the redshift of collapse, a constant (18-7T 2 ) for 
the f2 = 1 case and a otherwise a weak function of z 
fitted in Kitayama & Suto (1996). Navarro, Frenk, & 
White (1995) have shown that this relation with (3 = 
1.07 ± 0.05 is an accurate approximation to the results of 
N-body/hydrodynamic simulations in the case of an SI = 1 
universe. Here we follow Eke, Cole & Frenk (1996) and fix 
(3=1 for all cosmological models. 

Once the gas has fallen into a potential well, we adopt 
a simple model to calculate the time scale at which the 
gas will cool and form stars. Neglecting the gravitational 
energy associated with further gas infall we have 



n»(r)A(T(r),Z(r)) 



3 dT(r) 

—n (r )Kb 

2 K ' dt 



(11) 



where n e (r) is the number density of electrons at a radius 
f, n to t(r) is the total number density, T(r) is the tempera- 
ture of the gas, ks is the Boltzmann constant, and A is the 
radiative cooling function, which is strongly dependent on 
the metallicity Z{r) of the gas. Here we take A as tab- 
ulated by Sutherland and Dopita (1993) for equilibrium 
configurations, although strictly speaking these estimates 
are not exact due to the difference in abundance ratios 
for alpha elements in SNell outflows as compared to solar 
proportions. Considering only helium and hydrogen the 
mean molecular weight is pm p = p g /n = m p 4/(8 — 5Y) 
and the ratio of the electron and total number density are 

related by n = = |fe r~rrrs i where Y is the helium 

J 1 "tot 8(y 1 — 1)+3' 

fraction by mass, here taken to be 0.25. 



For a general configuration of gas, Eq. (11) is difficult 
to solve, and we therefore adopt a simple model for each 



collapsing halo as an isothermal sphere (p(r) 



oc r 



l ) with 



constant metallicity. We then follow the simple heuristic 
model of White & Frenk (1991) in which all the gas within 
some "cooling radius" r coo i cools instantaneously and all 
the gas outside this radius stays at the virial temperature 
of the halo, with r coo i moving outward with time. For an 
isothermal sphere this gives 



dM, 



cool 



dt 



dTc. 



dt 



= W^t 2 



A(T vir ,Z) 

10 — 2ll Cl'SS S~ 1 ( 



•^ 1/2 M yrTi2) 



where /hot is the fraction of the halo mass in the form of 
hot gas. While derived from a simple model, this expres- 
sion is in good agreement with similarity solutions for cool- 
ing flows given by Bertschinger (1989) (which are ~ 28% 
smaller than this expression), and the one-dimensional 
simulations of Forcada-Miro & White (1996) (which are 
~ 15% smaller). We therefore follow Somerville (1997) 
in multiplying the right hand side of this equation by an 
overall factor /o which wc take to be 0.8. Finally, we ap- 
proximate /hot — Afbaryon/Mhaio = O 6 /f2 at all times, 
and take the time at which M coo i = Mb ary0 n to be time at 
which a galaxy is formed. 

We show this formation time as a function of mass for 
four different metallicities at z = 1 in Figure |]. From 
this diagram we see that while cooling times are negligible 
for dwarf galaxies, cooling times for primordial-abundance 
halos with masses above ~ 1O 12 M are sufficiently long 
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to prevent galaxy formation at low redshifts. This is de- 
pendent on metallicity, as we saw in Figure |l|. In our 
simulation we therefore neglect cooling times for halos 
with masses < 1O 1O M0, while using the appropriate pre- 
enriched cooling times for more massive objects. 

Note also the sharp rise in cooling times at temperatures 
below w 10, 000 K, due to the lack of atomic transitions 
with energies below a few eV. Cooling at lower temper- 
atures can only occur through molecular transitions, and 
thus the collapse of T v ; r < 10 4 primordial gas clouds de- 
pends strongly on the existence of molecular hydrogen. A 
detailed investigation of the production and dissociation of 
H2 in the early universe has been conducted by Haiman, 
Rees, & Loeb (1997). As H 2 is easily photo-dissociated 
by 11.2-13.6 eV photons, to which the universe is other- 
wise transparent even before reionization, they find that 
a UV flux of <J 10~ 22 erg cm _2 s~ 1 Hz _1 sr _1 , is capable of 
dissociating all H2 in collapsing halos (see also Ciardi et 
al. 1999). As this is more than two orders of magnitudes 
smaller that the reionizing flux, we assume in our simu- 
lations that a small population of early stars quickly de- 
pleted the primordial gas of molecular hydrogen. Thus 
we exclude the formation of all galaxies with virial tem- 
peratures below 10, 000 K, which sets the lower halo mass 
value to « 5 x 10 7 M Q /i -1 , justifying the choice of spatial 
resolution for our simulations. 

3.4. Heating and Enrichment of the IGM 

Once a dwarf galaxy has been formed, we assume that 
a galactic outflow is immediately generated due to SNell 
from the short-lived massive stellar population. In order 
to get a rough estimate of the mass of metals injected into 
the ICM, the radius over which they are mixed into the 
surrounding gas, and the heating of the ICM, we model 
them according to the simple analytical model developed 
in Tcgmark, Silk, & Evrard (1993, hereafter TSE). Here 
the authors treat the outflow as an expanding shell that 
sweeps up most of the baryonic IGM and loses only a small 
fraction f m = 0.1 to the interior, such that the mass of the 
shell is 

m(t) = ^R{tfp b {\ - f m ) (13) 

where R(t) is the radius of the outflow. Acceleration of the 
shell is due to the internal pressure and deceleration from 
gravitational breaking, both estimated in the thin shell 
approximation (Ostriker & McKee 1988). With these ap- 
proximations R becomes 

5_ 3 • 2 H*R GM 

R -nJPR~R {R ~ HR) (14) 

Note that the gravitational term here includes a correction 
for the presence of a halo of mass M out of which the out- 
flow expands. The pressure of the hot gas within the shell 
is given by p = where Et is the thermal energy of the 
interior, and energy conservation gives = L — p^jr, 
where L is the luminosity, incorporating all heating and 
cooling of the interior plasma. 

The evolution of the blast wave is completely determined 
by the assumed luminosity history of the blast wave. Fol- 
lowing TSE we consider five contributions to this luminos- 
ity: the cooling by Compton drag against the microwave 



background (I com p), the cooling due to bremsstrahlung 
and other two body interactions (L ne 2) as in Eq. (|Tl|) , en- 
ergy injection from supernovae (L sn ), cooling by ionization 
of the IGM (Li on ), and heating from collisions between the 
shell and the IGM (Ldi SS ); such that, 

L = —L comp — L nc 2 + L sn — Lj on + Ldiss- (15) 

The first of these we fix according to Eq. (4) of TSE, while 
L ne 2 is negligible in this regime. The other parameters re- 
quire some calibration to the observed properties of out- 
flows from dwarf ellipticals and our overall model of the 
IGM at the time of dwarf formation. Let us consider these 
individually. 

Following TSE we assume that the energy in the blast 
wave is proportional to the number of SNell and winds 
for O and B stars and that this in turn is proportional 
to the baryonic mass which has collapsed. Energy in- 
jection occurs over a relatively short time during which 
the dwarf galaxy is undergoing a starburst phase. We as- 
sume that a fraction of €sf of the baryons in each galaxy 
collapse to form stars, that one supernova is formed for 
every 100 M© of stars in the halo, and that each su- 
pernova provides a typical energy output of 10 51 ergs, 
with an equal contribution coming from winds from mas- 
sive stars. If a fraction e W i n d of this input goes into the 
galactic outflow then the total energy from the super- 
novae is E = e W ind£sf2.0 x 10 55 ^-Af 6 - Assuming that this 
energy is released at a constant rate during a period of 
iburn = 1.0 x 10 8 years, we find 

L s „ = e wind e 8f 1.6 x 10 e L Q /hM 6 x -± (16) 

"0 

where the factor of h is a result of the units we have cho- 
sen for Mq. Dekel and Silk (1986) have proposed that the 
observed differences in surface brightness and metallicity 
between the observed classes of diffuse dwarfs and normal 
galaxies can be understood in terms of outflows occurring 
in halos with virial velocities below critical value on of the 
order of 100 km/s, which corresponds to approximately 
2.5 x 10 9 esF(l + z) 3 / 2 f2b/f2o supernovae. Here we take a 
conservative approach and exclude all bursts greater than 
2 x 10 7 esF^b/^0i which is typically a few hundred thou- 
sand supernovae and corresponds to a maximum luminos- 
ity of 

L sn , max = £windesf3.2 X 10 L© X — (17) 

"0 

achieved by objects above 2 x 10 9 M©. 

Another potential source of drag on the expanding shell 
is energy losses to ionizing the IGM, given by Li on = 
/neutral fm rib E AirR 2 [R — H R] , where rib is the number 
density of baryons, Eq ~ 13.6 eV, and /i on is the ioniza- 
tion fraction of the IGM into which the shell is expanding. 
As they were trying to develop a model for IGM reion- 
ization by winds, TSE assume /neutral = 1 at all times. 
Here we choose for our fiducial model the case in which 
blast waves expand into a medium that has completely 
been pre-ionized (/neutral = 0) possibly by ionizing pho- 
tons from active galactic nuclei, or from ionization fronts 
from the massive stars that precede the winds. 

Finally Ldi SS accounts for the heating of the interior 
plasma due to collisions between the shell and the IGM. 
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Taking this to be some fixed fraction /di ss of the kinetic 
energy lost by the shell gives 

^diss = /dissTTTT^ — HR) 3 . (18) 
lit 

In TSE the authors consider the extreme cases in which 
/diss = and /diss = 1) showing that their results are rel- 
atively insensitive to this parameter. Here we simply fix 
/diss = 0.5 throughout. 

With these approximations we can use the solution given 
in TSE to calculate the full expansion history of the galac- 
tic outflow, with their L* and appropriately modified to 
account for the gravitational potential of the collapsed halo 
as given by our Eq. (|l6|). We step forward in time (back- 
wards in redshift) in increments of Az — 0.01 starting at 
a maximum redshift of z max = 25, using a fourth-order 
Runge-Kutta technique. Figure || shows the first galaxy 
to undergo an outflow, a 1.1 x 10 8 M Q dwarf expanding 
into an ionized medium at z — 25 in one such simulation. 

As our simulations make no attempt to include ioniza- 
tion and do not account for the presence of ionizing ra- 
diation from stars within the dwarf galaxy or an external 
source of UV radiation, we expect our solution to only be 
accurate while the temperature of the material is greater 
than 10, 000 K (k b T = 13.6 eV). After this point, we as- 
sume as a reasonable approximation that the halo simply 
expands with the Hubble flow and remains at a fixed tem- 
perature of 10,000 K. 

In this figure, we also include for comparison the results 
of one such blast wave expanding into a neutral medium. 
This closely traces the ionized solution up until the tem- 
perature of the outflow drops to ~ 10, 000 K. 

Finally, we model the metallicity of the expanding halo 
by assuming that each supernova typically ejects 2Mq of 
metals into the intergalactic medium. This number is con- 
sistent with the average stellar yields in SNell simulations 
as compiled in Nagataki & Sato (1997), although there 
are significant theoretical uncertainties between the vari- 
ous simulations. The total mass of metal injected into the 
interstellar medium of the star forming galaxies is then 

M M = 0.02erf(n 6 /fi )M 6 . (19) 

Our model has two important free parameters, e s f and 
e-wind- While the amount of energy in the winds is directly 
dependent on the product of these two parameters, their 
ratio is somewhat more unconstrained. Our approach here 
then will be to fix the star formation efficiency at a some- 
what standard value of e s f = 0.1 and allow the fraction of 
the SN energy channeled into the wind to vary between 5 
and 20 percent to quantify model uncertainties. 

With this simple model we can construct from X Mi 
and the full histories of the blast waves expand- 

ing into the intergalactic medium: the center of each 
expanding bubble (X Wast = {x blast , xi? last , ...}), its co- 
moving radius as a function of redshift (lZ hlast (z) = 
{i? blast (z), i?2 last ( z ); ••■}), temperature as a function of 
redshift (T blast (z) = {"7\ blast (z), T, blast (z) , ...}), and mass 
of ejected metals (7W|' ast = {Af bfast , M|] ast , ...}). While 
lZ blast (z) and 7" blast ( z ) are updated in redshifts intervals 
of Az = 0.01, the history is only stored in intervals of 
Az = 0.05 to save memory. Note that the jagged ap- 
pearance at late times of the comoving size of the out- 
flow in Figure is caused by this discrepancy in times 



steps, rather than by a numerical instability. Note also 
that as outflows only occur in objects with masses <J 10 10 
for which cooling times are negligible, we can compute the 
entire outflow history of our Monte Carlo volume with- 
out considering the impact of metal enrichment on the gas 
cooling times. Suppression of galaxy formation due to out- 
flows from nearby objects, however, must be considered for 
outflow-scale second generation objects. 

3.5. Formation of 'Second Generation' Objects 

As we saw in §2, galactic outflows affect subsequent 
galaxy formation both by shocking the IGM as well as en- 
riching it with metals. Shocks can suppress the formation 
of nearby galaxies by two main mechanisms, heating the 
halo gas or striping it from the dark matter. In the first 
case the gas associated with a collapsing halo is heated to 
above the viral temperature of the halo but remains bound 
to the collapsing dark matter. The thermal pressure of the 
gas then overcomes the dark matter potential and the gas 
expands out of the halo, preventing galaxy formation. In 
the second case suppression is caused by stripping of the 
baryons by a shock from a nearby source. The momen- 
tum of the moving shock is sufficient to carry with it the 
gas associated with the halo of an unvirialized perturba- 
tion, thus emptying the dark matter halo of its gas and 
preventing a galaxy from forming. 

In a companion letter (Scannapieco, Ferrara, & Broad- 
hurst 2000), we evaluate both these scenarios of halo sup- 
pression and find the stripping mechanism to be domi- 
nant, because most shock- heated clouds can radiatively 
cool within a sound crossing time. As in the case of out- 
flow generation, momentum transfer between the outflows 
and the collapsing halos is a complicated process, depen- 
dent on the overall density profile of the collapsing halo 
and the possibility of smaller collapsed objects within it. 
Again, our approach in this exploratory study will be to 
adopt simple criteria to estimate on average when the mo- 
mentum in a shell is sufficient to remove material from a 
collapsing halo. 

We allow for baryonic stripping to occur only if a halo 
has not yet virialized, and assume that stripping occurs in 
all cases in which a shock moves through the center of the 
pre- virialized halo with sufficient momentum to accelerate 
the gas to the escape velocity. Thus we exclude all objects 
for which 

fM s v s >M c v e , (20) 

where / = £ 2 /4r 2 is the solid angle of the shell that is 
subtended by the collapsing halo, i is the radius of the 
collapsing region when the shock moves through its cen- 
ter, M s = 4^i7f,i7(z)/r2o/°c^s is the mass of the material 
swept up by the shock and M c is the baryonic mass of the 
collapsing halo. We can estimate the comoving radius of 
the collapsing halo as 

£ = i? Mj (l + <5 NL )" 1/3 , (21) 

where <5nl is the nonlinear overdensity of the region, which 
is given in the spherical collapse model as 

9 (0- sing) 2 
1 + * NL= 2 (l-coBg)3 ' (22) 
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where the collapse parameter 9 is given by (8 — 
sin(9) 2 / 3 7r~ 2 / 3 = D(z cmss )/D(z c ) where z cross is the red- 
shift at which the shock moves through the center of the 
halo. Note that in the elliptical collapse case, z c is com- 
puted according to Eq. (ji|), slowing down the condensation 
of the cloud somewhat. Note also that we treat halos with 
subcondensations that have already virialized in the same 
way. In this case the shock would lose energy imping- 
ing on the condensed sub-objects but would also have less 
uncondensed gas to accelerate to escape velocity. We as- 
sume for our purposes here that these effects roughly com- 
pensate, although more detailed hydrodynamical studies 
would help to refine this comparison. 

Finally, we calculate the metallicity of each collapsing 
sphere by volume averaging the contributions of each of 
the expanding blast waves that pass within the collapse 
radius r co \. Each halo is assigned a collapse mass in met- 
als M^ 1 kl which is taken to be zero initialy and modified 
by each blast front that passes within the collapse radius. 
For each such occurrence, the mass in metals is updated 
to 

M%i — M# + ^g*M|^ (23) 

3 n Mi 

where Overlap is the volume of intersection of the two 
spheres. By dividing this mass by the total baryonic mass 
of the galaxy we can compute the initial metallicity of the 
object. This value is then used to compute the collapse 
time of the object by using the appropriate cooling func- 
tion, A(T, Z). 

4. RESULTS 

In this section we summarize the results of our simu- 
lations, both in relation to the suppression of low-mass 
galaxy formation by baryonic stripping of pre-virialized 
halos, as well as the pre-enrichment of high mass galax- 
ies by metals carried in galactic outflows. We show that 
the tendency for small galaxies to disrupt the formation 
of their neighbors helps to explain the small number of 
Milky Way satellites relative to Cold Dark Matter model 
predictions, and leads to a lack of small-scale galaxy merg- 
ers, suggesting a bell-shaped luminosity function for ellip- 
tical galaxies. Pre-enrichment helps to explain the high 
G-dwarf metallicities found in our own and other galax- 
ies, but the shorter cooling times of metal enriched clouds 
seems to have little effect on the overall number of large 
galaxies formed. We examine the effect of varying model 
parameters on these features and discuss the compatibility 
between the star formations rates and IGM heating in our 
models and constraints from optical observations and mea- 
surements of thermal distortions in the cosmic microwave 
background. 

4.1. Shocks and Low-Mass Galaxy Formation 

Having constructed a distribution of halos with collapse 
times and spatial orientation in good agreement with ana- 
lytical expressions, we now consider the impact of outflows 
and enrichment on galaxy formation within this distribu- 
tion. We first restrict our attention to the ACDM model 
of Figure 1, as the currently favored cosmological model, 
and consider the effects of varying the cosmology and other 
model parameters in §4.4. In Figure || we show the col- 
lapse and formation of outflows in our fiducial model with 



Cwind = 0.1. Here we see that regions effected by galac- 
tic outflows are highly correlated with regions in which 
new objects are collapsing, illustrating the necessity of an 
inhomogencous approach. 

In Figure ^ we show the volume fraction in the expand- 
ing shells as a function of z. This quantity is not esti- 
mated on a cell-by-cell basis, but rather approximated by 
its value in the case in which the positions of the bubbles 
are uncorrelated, 

(47T y-ijV r>\Aast I y \&\ 
- 3 r j , J ■ ( 24 ) 
(12Mpc//i) 3 J V ' 

Note that this quantity overestimates the volume within 
outflows as the correlations between outflow positions can 
be significant, and are even greater in a full N-body ap- 
proach than in our simple Monte Carlo. 

Also shown on this plot is the volume-averaged temper- 
ature within the shells and the overall volume-averaged 
temperature of the simulation. Here we see that at red- 
shifts much greater than ~ 5, the areas near dwarf star- 
burst galaxies were significantly hotter than the IGM as 
a whole. This discrepancy is somewhat reduced at later 
times, when the outflows begin to fill a significant frac- 
tion of the total volume of the universe. Even at these 
late times, however, the universe remains inhomogeneous, 
with new starbursting dwarf galaxies forming within re- 
gions only mildly impacted by the earliest outflowing ob- 
jects. Thus the root mean squared (RMS) fluctuation in 
the temperature remains significant at redshifts p> 1. 

This high degree of inhomogeneous shocking naturally 
has a large impact on galaxy formation at masses below 
~ 10 10 M Q . In Figure || we plot the number of collapsed 
halos as a function of redshift, along with the number of 
unsuppressed objects as described in §3.4. While the earli- 
est dwarf galaxies form along with the collapsing halos, by 
z » 12 outflows become important. As M ^ 10 9 M Q halos 
collapse after this redshift, galaxy formation at this scale 
is highly suppressed, leading to a "mass desert" between 
the major population of polluting dwarf galaxies and the 
larger galaxies that are seen today. It is also clear from 
these figures that galaxy formation is likely to have oc- 
curred in two stages, with the dwarf outflow population 
forming mostly at the highest redshifts, and larger galax- 
ies forming only later when halos of sufficient mass began 
to collapse. 

Also in this figure, we plot the total number density 
of collapsed halos with masses between 1.1 x 10 8 M Q and 
1.4 x 10 11 M Q , along with the total number of galaxies in 
this mass range. Notice that only ~ 30% of the halos are 
populated, consistent with the factor of ~ 4 suppression 
needed to reconcile the number of predicted and observed 
Milky- Way satellites (Kauffmann, White, & Guiderdoni 
1993; Klypin et al. 1999; Moore et al. 1999). Our scenario 
thus provides a natural mechanism for the formation of 
"dark halo" satellites around our galaxy, which may be 
associated with the abundant High- Velocity Clouds as dis- 
cussed by Blitz et al. (1999). 

4.2. Enrichment and High-Mass Galaxy Formation 

In Figure ^ we plot the volume-averaged metallicity 
within the outflows and the simulation overall, as well as 
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the RMS variation between the shells. Here we see that 
like the temperature, the mean metallicity within the out- 
flows is the greatest at early times, and decreases to a few 
times the mean metallicity of the universe at a redshift of 
~ 5. The mean IGM metallicity from z = 1 to z = 5 ranges 
from about log 10 [Z/Z Q ] ~ —1.8 to log 10 [Z/Z Q ] ~ —2.6 to 
in good agreement with observed metallicities in Ly ab- 
sorption systems (e.g. Cowie & Songaila 1998, McDonald 
et al 1999, Ellison et al. 1999), as well as previous studies 
of enrichment by dwarfs (Nath & Trentham 1997). The 
RMS variations in the metallicities within the outflows are 
even more severe than the temperature variations however, 
and remain several times larger than the mean metallicity 
at all times. 

This highly inhomogeneous distribution of metals is con- 
sistent with observed metallicity inhomogeneities as dis- 
cussed in §2.2 and suggests that the mean metallicity of 
the gas out of which most galaxies formed may be signif- 
icantly higher than that of the IGM. Figure || also shows 
the initial metallicity of the IGM newly formed galaxies. 
Note that this does not include the contribution from the 
metals in the progenitor galaxies themselves. Here we see 
that indeed the mean metallicity of collapsed objects is 
an increasing function of mass which is expected to be 
log 10 [Z/Z Q ] ~ —1 for Milky- Way sized objects, and higher 
than the mean IGM metallicity at for all mass scales. This 
value is comparable to the ~ O.IZq pre-enrichment nec- 
essary to explain the lack of low-metallicity G-dwarfs in 
the Solar neighborhood. The overall relatively high initial 
metallicity is also consistent with observations that sug- 
gest the need for a "floor" of around O.IZq in the metal- 
licity distribution of stars in other massive galaxies. Note 
that we make no attempt to trace subsequent star for- 
mation and interstellar medium (ISM) enrichment, and 
our values only provide a lower limit on the lowest metal- 
licity stars in disk galaxies that formed largely from the 
IGM gas. Note also that these metallicities are slowly in- 
creasing with mass, suggesting an alternative explanation 
for the mean metallicity luminosity relationship, which is 
usually understood as ISM enrichment within galaxies and 
their progenitors (see eg. Vader 1986; Ferguson & Binggeli 
1994; Kauffmann & Chariot 1998). 

As collapsing halos in which galaxy formation is sup- 
pressed are found closest to outflowing galaxies, where 
both shock velocities and metallicities are high, it is likely 
that whatever residual gas remains in these objects would 
be of similar of even higher metallicity than that of un- 
supressed objects. While the calculation of these metal- 
licities is beyond the scope of our simulations, it is inter- 
esting to note that some of the high velocity clouds may 
also have metallicities that actually exceed those of many 
dwarf galaxies although may still be somewhat less than 
suprasolar values predicted by a "galactic-fountain" model 
(Sembach et. al 1999; Wakker et al. 1999). 

In Figure [l(] we show the effect of cooling on galaxy for- 
mation in large-mass halos. Here we see that while halo 
collapse and galaxy formation are almost simultaneous at 
masses ^ 5 x 10 10 M Q , long cooling times suppress all 
galaxy formation on scales > 5 x 10 12 M Q . The increased 
initial metallicities help to accelerate galaxy formation in 
objects between these mass limits somewhat, although the 
overall final number densities remain unchanged. 



4.3. Properties of Elliptical Galaxies 

As elliptical galaxies tend to be found in the most dense 
and enriched regions of space, it is natural to expect that 
outflows would have had the largest impact on these ob- 
jects. In order to study this connection, we follow the con- 
ventional wisdom that ellipticals correspond to mergers of 
(see e.g., Barnes 1992; Hernquist 1993) large progenitors. 



We therefore define M, 



Mi 



as the total contribution 



merger, k 

to the mass of halo k of mass Mi due to mergers of unsu- 
pressed galaxies of mass scales Mj_i and M,_ 2 , that is 



M. 



M, 



M(_i 


A/z-2 





if 



if 



otherwise 

r Mi 



Mi. 



otherwise 



Ms_ 2 



< Rm, 



(25) 



where we consider only objects that have cooled and have 
not been swept away by shocks while forming. We then 
identify as ellipticals all collapsed objects with M m ^ k > 
0.5Mi, in which half of the total mass comes from a merger 
of large progenitors. 

We fix this threshold in order to approximate the ob- 
served field elliptical fraction of ~ 15% (Baugh, Cole, & 
Frenk 1996), and ask the question of what mass scales 
correspond to these objects. 

In Figure [ll] we plot the total space density of galaxies 
in our simulations, along with the space density of galaxies 
identified as ellipticals. In this figure we see that the mass 
distribution of ellipticals is quite different than that of the 
overall galaxy population. As the outflows from dwarf 
galaxies suppress the formation of nearby dwarfs, there 
are very few mergers at small masses, with virial temper- 
atures much smaller than that typical of outflows. Thus 
while the number density of halos continues to increase 
with decreasing mass, the number density of ellipticals is 
the same at 1.4 x 10 11 M Q and 4.4 x 10 10 M Q , and decreases 
dramatically for even lower masses. 

Using the observed Faber- Jackson relation between the 
luminosity of elliptical galaxies and the velocity disper- 
sion, we can recast this as an absolute r-band magnitude of 
M r = —18 — 10 log 10 ctioo, (Oegerle & Hoessel 1991) where 
(Tioo is the velocity dispersion in units of 100 km/s. This 
gives the luminosity function of elliptical galaxies shown in 
the bottom of Figure [ll]. While these results are necessar- 
ily crude, and are bounded at the high-luminosity end by 
the finite size of out simulations, they nevertheless natu- 
rally reproduce the observed bell-shaped luminosity func- 
tion of elliptical galaxies (see eg., Bromley et al. 1998). 

4.4. Effects of Model Uncertainties 

We have only two free parameters in our simplified mod- 
eling, both relating to the gas outflow, and all other model- 
ing based on a simple spherical or elliptical collapse model, 
involving the usual assumptions. Our conclusions are most 
strongly dependent on the product of the two free param- 
eters e w i n d and e s f , which together determine the total en- 
ergy being channeled into galactic outflows. In our fiducial 
model, we took e 8 f = 0.1 and e w i n( je s f = .01. These canoni- 
cal values are motivated by the current observational work 
on dwarf galaxy outflow energetics, but a large spread 
is naturally expected and thus we have some freedom in 
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changing this efficiency. We therefore examined two ex- 



treme cases with e v 



= .05 and e„ 



= .2 to asses the 



robustness of the features described above. These cases 
bound the shaded regions of the results plotted in Figures, 
7, 8, 11, 12. 

In Figure ^ we see the impact of model uncertainties 
on the total volume fraction and temperatures of outflows 
in our simulations. The impact of these differences on 
the suppression of dwarf galaxies is shown in the shaded 
regions in Figure |^. Here we see that the existence of two- 
stage evolution with a "mass desert" in the 10 9 — 10 10 M© 
range persists for all values of these parameters. The to- 
tal number of unsuppressed dwarf galaxies ranges from 
10% — 30% reinforcing the idea that outflow suppression 
may be important in reproducing the small number of 
Milky- Way satellites. 

The effect of varying e s f e w ind on the distribution of ellip- 
tical galaxies is shown in Figure 11. While the uncertain- 
ties in this case are somewhat more severe than in Figure 
||, a large fall-off in the number of objects with masses 
below ~5x 1O 1O M0 is present for all cases. 

We show the metallicities of the IGM component of the 
formin g ga laxies in each of these cases in the upper panel of 
Figure |2[ In both cases the mean metallicity of collapsed 
objects is an increasing function of mass which is higher 
than the mean IGM metallicity. The predicted Milky- Way 
scale IGM metallicity is again always log 10 [Z/Z Q ] w —1 
and thus the G-dwarf problem is easily understood. 

In the lower panel of this figure, we compute the star 
formation rate in dwarf galaxies to assuage any concern 
that our code might be pumping an inordinate amount 
of energy into the IGM. Note that our models do not in- 
clude subsequent star-formation in disk galaxies and thus 
this figure represents a lower bound to the number of stars 
formed at each redshift. Nevertheless, our values are well 
within observational constraints (see eg., Adelberger & 
Steidel 2000). 

We can also compare our models to the COBE con- 
straints on the optical depth to the surface of last scatter 
as well as the the distortion of the cosmic microwave back- 
ground (CMB) spectrum due to the Sunyaev-Zel'dovich 
effect. We can estimate the optical depth in our model as 

r = 5.9 x HT 3 ^ f ^ dz^n(z) (1 + z) 2 F(z), (26) 
ll J dz 

where x is again the comoving distance defined such that 
aoHo = c. This analysis gives r = 0.05 for the fiducial 
model, and varies from r = 0.04 to r = 0.06 in the low 
and high energy outflow cases, which are all well within 
the observational limit of r <; 0.5 (Griffiths, Barbosa, & 
Liddle 1999). 

The degree of CMB spectral distortions is given by the 
Compton-y parameter which is the convolution of the op- 
tical depth with the electron temperature along the line of 
sight (Zel'dovich & Sunyaev 1969; Sunyaev & Zel'dovich 
1972). Thus 

y = 1.0 x 10" 8 -^ / dz-fT 5 (z)n(z) (1 + z) 2 F(z), 

where T$(z) is the mean temperature within the outflows 
as a function of redshift in units of 10 5 K. In this case 



y is 1.6 x 10~ 6 in the fiducial model, and varies from 
1.0 x 10~ 6 to 2.3 x 10~ 6 for the extreme cases. These 
values are again within the observational constraint of 
y < 1.5 x 10 -5 (Fixsen et al. 1996). Nevertheless, our mod- 
eling suggests that the spectral distortions due to galactic 
outflows should be observable with the next generation of 
cosmic microwave background experiments. 

Varying the second parameter in our model, e s f, while 
keeping e s fewind fixed acts as an overall shift in the star 
formation rate and metallicities in Figures | and |l| and 
can be estimated simply by eye. Including the ionization 
drag term in Eq. (11) and allowing the shocks to act as the 
ionizing source of the IGM has little effect on our results 
as suggested by Figure || and verified by a full simulation. 
Thus we see that all the major features of the fiducial 
model persist over a wide range of parameters. 

We have also examined the cosmological dependence 
of our results by examining a flat model with parame- 
ters as described in §2.3 and e w i n d = 0.1. In this case 
only ~ 50% suppression of galaxies occurs, again with the 
most impact on galaxies in the few times 10 9 Mq range. In 
this model elliptical formation falls off below the 1O 11 M 
scale but more gradually than in the open case, and initial 
Milky- Way metallicities again about log 10 [Z/Z@] = —1.0. 
Thus even in a model in which structure forms much more 
quickly than suggested by observations (see eg., Bahcall, 
Fan, & Cen 1997), the major features in our model persist. 

5. DISCUSSION 

Our treatment of galactic outflows is intended as an ex- 
ploratory study of these processes, and makes a number 
of approximations that should be made explicit. While a 
Monte Carlo approach uncovers many of the issues of inho- 
mogeneity that can not properly be studied analytically, 
it fails to capture the nonlinear structure and clustering 
present in a full N-body treatment. At late times peaks 
will be grouped along sheets and filaments, and thus clus- 
tering is likely to be more severe than in our simulations 
and of a more complicated nature. 

A number of uncertainties also arise from our simple 
treatment of IGM enrichment and shocking. One concern 
is the presence of radiative feedback from ionization fronts 
around dwarf galaxies. Our assumption that fronts pre- 
cede the formation of outflows, and thus each shell can 
be treated as expanding into an ionized medium is prob- 
ably a reasonable one, while our placement of a "tem- 
perature floor" at the ionization temperature of hydrogen 
at all times in our simulations is stronger approximation. 
Fortunately the extreme fragility of molecular hydrogen 
suggests that the formation of objects with virial tempera- 
tures below ~10,000 K may have not been largely affected 
by reionization. While the inhomogeneous structure of 
reionization is likely to have had a great impact on the 
CMB fluctuations (Aghanim et al. 1996; Miralda-Escude, 
Haehnelt, & Rees 2000; Scannapieco 2000), the formation 
of 10 7 M e galaxies is likely to have been halted by the 
dissociation of molecular hydrogen by the first stars, long 
before reionization took place. 

A bigger concern may be the structures of outflows and 
the ejection fractions of interstellar gas and metals. This 
has been studied in detail by Mac Low & Ferrara (1999) 
and Ferrara & Tolstoy (2000), who conclude that efficient 
ejection of the interstellar medium or "blowaway" occurs 
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only in halos with masses <J 1O 7 M0. The issue of whether 
dwarfs retain a sizeable fraction of their gas, however, 
is to some degree decoupled from the formation of out- 
flows. In halos in the mass range 10 7 Af Q <; M <, 10 9 M Q , 
a "blowout" occurs in which the super-bubbles around 
groups of SNell punch out of the galaxy, shocking the 
surrounding IGM and efficiently ejecting metals while fail- 
ing to excavate the interstellar medium of the galaxy as a 
whole. Such a scheme in which outflows are formed pri- 
marily of the enriched IGM surrounding dwarf galaxies is 
equivalent to outflows of galactic gas for the purposes of 
our simulations. Additionally such a picture would help to 
reconcile observations of expanding high metallicity shells 
around dwarf galaxies as described in §2 with observa- 
tions of multiple episodes of and HI gas in dwarf spheroidal 
galaxies (Smecker-Hane et al. 1994; Grebel 1998) some of 
which even suggest that many of these objects are gas-rich, 
but with extended HI envelopes (Blitz & Robishaw 2000). 

That being said, the TSE model adopted for outflow 
evolution in our simulations is almost certainly oversim- 
plified. The "blowout" scenario described by Mac Low 
& Ferrara (1999) ejects matter perpendicular to the thick 
disk, and is thus quite different from the spherical out- 
flow modeled here. Even if the winds from dwarf galaxies 
can be reasonably approximated as spherical shells at early 
times, such outflows of low-density heated gas would neces- 
sarily become Ray leigh- Taylor unstable as they expanded 
into the denser IGM. Thus material is most likely to flow 
out in a number of heated clumps at large distances. 

Also ignored in our model are shell-shell interactions 
which would contribute additional heating but slow ex- 
pansion in regions which two shocks meet. Similarly, an 
additional source of IGM enrichment is the merger mech- 
anism described in Gnedin & Ostriker (1997) and Gnedin 
(1998), in which a significant fraction of the interstellar 
medium of merging galaxies is ejected. 

Finally, the criteria for halo suppression used in our sim- 
ulation is over-simplified. While regions of space in our 
simulations that meet the stripping criterion, Eq. (p0[), are 
likely to have delayed or suppressed galaxy formation, the 
one-to-one correspondence assumed in our simulations is 
unlikely. In reality the formation of galaxies is likely to be 
a complicated function of the number and strength of the 
shocks moving through the regions, and the structure of 
these regions during and after halo collapse. This subject 
merits further investigation and would help to sharpen our 
conclusions. 

6. CONCLUSIONS 

While the impact of preheating and enrichment on the 
observed properties of galaxy clusters has long been rec- 
ognized, the impact on the observed properties of galaxies 
themselves has been little explored. By accounting for the 
observed outflows from dwarf galaxies, we have been able 
to show how many of the unexplained properties of galax- 
ies and the IGM can be naturally understood. 

Firstly, the suppression of low-mass galaxy formation by 
outflows provides a natural explanation for the factor of 
~ 4 discrepancy between the number of observed Milky- 
Way satellites and predictions from standard CDM models 
that do not include outflows. Suppression also provides a 
natural mechanism for the formation of "dark halos" which 
may be associated with the High- Velocity clouds. 



Secondly, baryonic stripping results in a bell-shaped lu- 
minosity function of ellipticals. The lower the mass of a 
halo, the more likely it is to generate an outflow that strips 
material from a similar mass neighboring pre-virializcd 
halo that would otherwise later form into a galaxy. This 
results in very few pairs of neighboring low-mass galaxies, 
and hence a relative deficit of major low mass-mergers. 

Finally, our models of enrichment of protogalactic gas 
predict a trend of increasing metallicity with galaxy mass 
in good agreement with inferences from observations. The 
initial metallicity predicted for a Milky- Way mass galaxy 
is ~ O.IZq providing a natural initial floor at the level re- 
quired to solve the G-dwarf problem and the more general 
lack of low metallicity stars in well studied massive ellip- 
tical galaxies relative to "closed-box" models of chemical 
enrichment. 

These results are persistent over a wide range of model 
parameters and cosmologies, and are not a result of fine- 
tuning parameters or invoking additional physics. 

While galaxy outflows are already incorporated into 
modern studies of galaxy formation, this is done only as an 
internal modification, ignoring pre-enrichment. Keeping 
track of the effect of outflows on neighboring halos is es- 
sential in understanding the properties of galaxy clusters, 
and hence it is not surprising that these effects would play 
a major role in the formation of galaxies as well. While the 
details await further investigation, it is clear that any com- 
plete picture of galaxy formation must account for heating 
and enrichment. 
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Fig. 1. — The impact of galactic outflows on galaxy formation. The solid lines show the lowest z at which a halo can collapse and still form 
a galaxy by z = 0, illustrating the enhanced cooling of metal-enriched virializcd objects versus objects formed of primordial gas. In each pair 
the upper line corresponds to primordial gas, and the lower line to gas that has been enriched to 0.1 Zq. The dotted lines show the collapse 
rcdshift of the 2cr (upper line) and la peaks (lower line), indicating the rcdshift at which objects of various mass scales form. The dashed 
line corresponds to a fixed virial temperature of 5 X 10 5 K, typical of outflows, and below which outflows are capable of suppressing galaxy 
formation. The upper panel lines are calculated in a J!o = 1 CDM universe (Qo = 1, = 0, fij, = 0.07, erg = 0.6, and T = 0.44) and the 
lower panel lines are for a ACDM model (Q = 0.35, tt A = 0.65, U b = 0.06, ct 8 = 0.87 and T = 0.18). 
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Fig. 2. — Check of our Monte Carlo space densities. The solid lines correspond to the 3-D Monte Carlo space density of halos at each 
redshift and the dotted lines correspond to the analytical Press-Schechter number density. The panels are labeled by mass in units of Mq. 
Here we see that over the reliable range of redshifts our method is in reasonable agreement with analytical predictions on all scales. 
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Fig. 3. — Comparison of analytical (solid lines) and Monte-Carlo (dotted lines) correlation functions at various mass scales labeled in units 
of Mjodot. The objects are calculated from the spherical collapse model and at redshifts such that v = 8 c Do/a(R)D(z) = 1.5. While both 
the analytical and numerical curves must be different at small values of r because of the manner in which they are calculated, good agreement 
exists at intermediate distances for all mass scales. 




Fig. 4. — Cooling times adopted in our calculations as a function of halo mass and virial temperature at a formation rcdshift of z = 1, in 
a flat cosmology. The solid line is for primordial abundances ([Fe/H] = -3), the dotted line corresponds to [Fe/H] = -2, the short-dashed line 
to [Fc/H] = -1, and the long-dashed line to solar abundances. 
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Fig. 5. — Outflow from a 1.1 X 1O 8 M0 galaxy with e w ; nt j = 0.1. The solid and dotted curves show the expansion into an ionized and neutral 
medium respectively. Note that the neutral fraction of the IGM is largely unimportant in determining the evolution of the outflow. 
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Fig. 6. — A slice though the simulation volume at four different redshifts. The filled circles show the radii of halos that have collapsed, 
while the open circles show the size of the galactic outflows that are in progress. Outflowing regions are highly correlated with regions in 
which new objects are collapsing at all redshifts, illustrating the necessity of a three-dimensional simulation. 
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Fig. 7. — Top: Volume filling factor, F(z) as a function of redshift. Bottom: Mean temperature within the outflowing shells (solid line) and 
the simulation overall (dotted line). The RMS fluctuations in the temperature of outflows is given by the dashed line. All lines correspond to 
the fiducial e w i nt j = 0.1 model, while the grey regions are bounded from above by the high-energy outflow case (e w i n d = 0.2) and from below 
by the low-energy outflow case (e wind = 0.05). 
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Fig. 8. — Suppression of low-mass galaxy formation by galactic outflows. In each panel, the upper line corresponds to the Monte Carlo 
space density of halos and the dashed line corresponds to the density of galaxies satisfying the criteria that the virial temperature be greater 
than the volume- averaged temperature of the collapsing halo, for the fiducial model. The shaded region is bounded from above by the 
low-energy outflow case and from below by the high-energy outflow case. The panels are labeled by mass in units of Mq. Outflows are most 
important in suppressing galaxy formation in the mass range 1O 9 M0 — 10 10 Mq. Note that we exclude the lowest mass scale from the total 
as it undergoes suppression when its viral temperature drops below 10, 000K at low redshifts, and the formation of these objects is therefore 
much more uncertain than at larger mass scales. 



18 








[SI 



&0 

o 




mass/M e 



Fig. 9. — Top: Mean mctallicity within the outflowing shells (solid line) and the simulation overall (dotted line). The RMS fluctuation in 
mctallicity within the outflows is given by the dashed line. Note that there is a large scatter in outflow metallicities, and the mean metallicity 
within the outflows is significantly higher the in the overall IGM average at all rcdshifts. Bottom: Mean metallicity of the IGM component 
of galaxies at the time of formation bounded by the RMS variations in galaxy metallicity. 




Fig. 10. — Suppression of high-mass galaxy formation due to long cooling times. The solid lines correspond to the Monte Carlo number 
density of halos at each redshift, the dashed lines correspond to the number density of cooled objects as calculated from primordial abundances, 
and the dotted lines to the number of cooled objects with metal enrichment. The panels are labeled by mass in units of Mq. Outflows accelerate 
galaxy formation on the ~ 10 12 Mq scale, however the overall number of galaxies at z = remain relatively unaffected. 
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Fig. 11. 



Top: Formation of elliptical galaxies. In each panel, the upper line corresponds to the Monte Carlo number density of halos at 



each redshift, and the lower line corresponds to the number density of galaxies with M n 



> 0.5M in the e w i n( j = 0.1 model. The shaded 



region is bounded from below by the e w i n d = 0.2 case and from above by the e w i n( j = 0.05 case as discussed in §5.4. Bottom: Elliptical 
luminosity function for our simulations. The solid line corresponds to the total r-band luminosity function for the elliptical galaxies in the 
fiducial model and the shaded region is bounded by the low and high energy outflow models. 
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Fig. 12. — Top: The mean metallicity of galaxies bounded by the RMS variations in galaxy mctallicity. The upper panel corresponds to 
the e w i n( j = 0.05 case and the lower panel corresponds to e w i n d = 0.2. Bottom: Star formation rate powering galactic outflows. Again the 
solid line is the fiducial model and the shaded region is bounded by the low and high energy outflow cases. 



